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Abstract 

We explore beyond Standard Model (bSM) physics signatures in the I + jets channel of if 
pair production process at the Tevatron and the LHC. We study the effects of bSM physics 
scenarios on the top quark polarization and on the kinematics of the decay leptons. To this 
end, we construct asymmetries using the lepton energy and angular distributions. Further, we 
find their correlations with the top polarization, net charge asymmetry and top forward back¬ 
ward asymmetry. We show that when used together, these observables can help discriminate 
effectively between SM and different bSM scenarios which can lead to varying degrees of top 
polarization at the Tevatron as well as the LHC. We use two types of coloured mediator mod¬ 
els to demonstrate the effectiveness of proposed observables, an s-channel axigluon and a 
M-channel diquark. 
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I. INTRODUCTION 


Most experimental observations at particle accelerators fit the Standard Model (SM) 
very well. However, there are some major puzzles to be solved. One needs to have 
physics beyond the standard model (bSM) to explain the presence of dark matter, ex¬ 
plain quantitatively the observed baryon asymmetry of the universe, and to explain the 
puzzle of dark energy. Looking for signs of bSM, one finds that most of the terrestrial 
experimental observations that are in tension with the SM results are in the properties 
of third generation fermions. For example B ^ tv [1] , h ^ fiT [2] and bb forward 
backward asymmetry (AFB) at LEP and Tevatron [3-5] show such deviations. One of 
these long standing puzzles is the top-quark AFB measured by the DO and CDF de¬ 
tectors at the Tevatron collider in 2008 [6, 7]. These observations by two independent 
collaborations were updated with full data from the Tevatron and were consistent with 
each other and in tension with the SM calculations until 2015. Recent experimental 
results from DO [3] and theoretical calculations [u] point towards the possibility that the 
anomalous nature of these observations may be a statistical phenomenon. 

Due to its large mass which is close to the electroweak scale and the implied con¬ 
nection with electroweak symmetry breaking, the top quark is an important laboratory 
for various bSM searches at colliders. In fact various proposals put forward to solve 
the different theoretical problems of the SM often involve modifications in the top sec¬ 
tor. Various extensions to the SM have also been proposed inspired by the possibly 
anomalous value of measured top-AFB at the Tevatron. These bSM proposals involve 
explanation of the AFB in terms of processes involving a) s-channel resonances like 
the axigluon, KK gluon, coloron [10-18] or b) t-channel exchange of particles with dif¬ 
ferent spin and SM charges like the Z' , diquarks etc [14, 16, 17, 19-25]. Effective 
operator approach has also been used in this context [26-28]. Measurements of other 
related observables such as lepton angular asymmetries and tt invariant-mass depen¬ 
dence of the top quark AFB are also compatible with the hypothesis of a heavy bSM 
particle, see for example [29, 30]. 

In this study, we will focus on the lepton -i- jets final state (pp/pp tt ^ bvlt) of the 
tt pair production process. This channel has a larger cross-section as compared to the 
dilepton -i- jets channel, and it has a much smaller background compared to the all jets 
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channel. For lighter quarks, hadronization smears the information available about their 
spin and polarization. The mass of the top quark is large enough that it decays into 
its daughter particles before strong interactions can initiate the hadronization process. 
Hence top-quark polarization leaves a memory in the kinematic distribution of the de¬ 
cay products and can be tracked [31,32]. We study the correlations between various 
kinematic asymmetries and polarization to distinguish between different sources of 
these asymmetries within an s-channel (axigluon) and a f-channel (diquark) extension 
of the SM. For the Tevatron, the top pair production process is dominated by qq colli¬ 
sions and at the LHC it is dominated by gg collisions which means that new physics 
can manifest at differently at the two colliders. 

A wide variety of observables have been studied in the literature to explore the 
top sector as a bSM portal [11, 16, 23, 33-37]. A brief review of some of these ob¬ 
servables which have been experimentally measured and are relevant to this work 
is presented in section II. In section III,IV we describe the flavour non-universal ax¬ 
igluon and diquark models which we use as templates for our analysis. Constraints 
on these models from top pair production cross-section and forward backward asym¬ 
metry at Tevatron, charge asymmetry, top quark pair production, dijet and four jet pro¬ 
duction cross-sections at LHC are discussed in section V. In section VI we construct 
the asymmetries which we use to explore the bSM models. In section VII we present 
the correlations between various asymmetries and discuss the role of top quark po¬ 
larization and kinematics in discerning the various regions of parameter space of the 
bSM models. We contrast our results for the axigluon and diquark models and the 
resulting conclusions can be generalized to other new physics scenarios. Our results 
are presented for the Tevatron ^/s = 1.96 TeV and the LHC ^/s = 7 TeV , 13 TeV We 
discuss the effects of transverse polarization coming from the off-diagonal terms in the 
top-quark density matrix in section VIII and then conclude in section IX. 

II. STATUS OF EXPERIMENTAL RESULTS 

We begin by summarizing some of the experimental results from the LHC and the 
Tevatron concerning the top quark and compare them with the corresponding SM cal¬ 
culations from the literature. 
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The measured tt production cross-section for the Tevatron at = 1.96 TeV is 

^Tmatron ^ 7 ^ g 4 ^ |■ 33 j |_|_|Q at = 7 TeV iS - = 

pp^tt t' l J V pp^tt 

173.30 ± 10.10 pb [39, 40]. These agree with the calculated SM NNLO cross-sections 
^ Tevatron, and = 177.30 ± 10.63 pb [42] for 

the LHC, within la. The uncertainties coming from the top-quark mass dependence of 
tt cross-section [43] have been included in the given LHC cross-sections. In the cal¬ 
culations in following sections, we use a common K factor for the bSM-i-SM to estimate 
the NNLO total cross-section. For the Tevatron, the K factor is Kjevatron = 1-39+q^q 
[44]. The K factor for the LHC is calculated using the NNLO cross-section cited above 
and LO cross-section calculated using CTE06I parton distribution functions (pdf) with 
factorization scale Q = 2mt. The errors in the K factors represent pdf uncertainties, 
scale dependence and statistical errors in the NNLO cross-section. For the LHC with 

v/^ = 7TeV,KL„C7 = 2.20+S;;4. 

The cross-sections impose a constraint on any new particle to have small couplings 
with the top quark and/or have a sufficiently high mass. It is interesting to note that we 
can still find a range of couplings of the bSM large enough to explain the measured 
anomalous top-quark and lepton asymmetries reported at the Tevatron and remain 
compatible with the measurements at the LHC. The AFB of the top quark in the tt 
centre-of-mass (CM) frame is defined as 

^ _NF{yt-yt>0)-NB{yt-yt<0) ,,, 

^FonvardBackward - > 0) + NB{yt " y? < 0) ^ ^ 

N{cos 9t > 0) — N{cos 9t < 0) 

N{cos9t > 0) + N{cos9t < 0)' 

where yt, yj are respectively the rapidities of the t and the t and 9t and 9i are their 
respective polar angles measured with respect to the beam direction. 

The CDF measurement of the tt CM frame AFB with the full data set is = 
0.164 ±0.045 [45]. The corresponding SM result is 0 at tree level in QCD. At NLO in 
QCD, the value predicted is 0.0589+gg^4g (the errors only represent scale variation) 
which upon including NLO electroweak corrections becomes, 0.0734+QQg5g [46]. Re¬ 
cently, AFB has been calculated at NNLO to be 0.0749+QQQgg in pure OCD and 0.095 ± 
0.007 including EW corrections [46] and including effective N^LO OCD, = 

0.100 ± 0.006 [u]. DO has come out recently with a measurement Afb = 0.106 ± 0.03 
[p] which agree with the theoretical results. However for the purpose of this study. 
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we use the CDF measurement which is still in tension with the SM and with the DO 
measurement. 

Since the LHC is a pp collider, its symmetric initial state makes the forward and 
backward regions trivially symmetric. For the LHC, instead of top quark AFB, a charge 
asymmetry (AC) is defined in the lab frame as. 


N(A|yf| 

V 

1 

o' 

V 

N(A|y,| 

V 

+ 

o' 

V 


(3) 


where A|i/t| = \yt\ - \yt\- The AC at the LHC is much smaller than the AFB at the 
Tevatron both in the case of the SM and of the bSM models aimed at explaining the 
Tevatron’s anomalous AFB. The measured value of AC with CMS and ATLAS combi¬ 
nation is Ac = 0.005 ± 0.009 [47]. The theoretical results for the SM values of the AC 
[48] (QED-i-EW-i-NLO QCD) are given in table I for different energies at the LHC. 


v^(TeV) 

Ac 

7 

0.0115(6) 

12 

0.0068(3) 

13 (From-Fit See Appendix A) 

0.0063 

14 

0.0059(3) 


Table I: Charge asymmetry in the lab frame at the LHC, as defined in eqn (3). 


Measurements have also been made for a number of other observables including 
Mff, rapidity-dependent top AFB [45], lepton and di-lepton asymmetries [49, 50], some 
of which show a deviation from the standard model [35] of up-to ^-3 a . Some CDF 
results are shown in table II and DO results [50] in table III. 

tt spin correlations have been measured using decay particle double distributions 
in polar and azimuthal angles at the Tevatron [51, 52] and the LHC [53, 54]. The po¬ 
larization of the top quark, as defined in eqn (12), has also been observed at CMS, 
for the LHC 7 TeV run to be 0.01 ± 0.04 [54] compared to the corresponding SM pre¬ 
diction from MC@NLO [55] 0.000 ± 0.002. The ATLAS collaboration also observed the 
polarization at 7 TeV beam energy, assuming CP conserving tt production and decay 
process, to be 0.035 ± 0.040 [56], in agreement with the SM prediction. 
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Asymmetry 

Experimental Value 

SM calculation 

Apg {or Aq^) 

U.U5'U_Qg26 

0.038 ± 0.003 

.Mtt>450GeV 
^ft FB 

0.295 ± 0.058 ± 0.031 

0.100 ± 0.030 

^Mti<450Gt;V 
^ft FB 

0.084 ± 0.046 ± 0.026 

0.047 ± 0.014 

aAI- 

^FB 

0.094 ± 0.024+°;°22 

0.036 ± 0.002 


Table II: CDF lepton and Mf^ dependent top level asymmetries [35, 45, 49] 


Asymmetry 

Experimental Value 

SM calculation 

Apg (or AeJ (extrapolated) 

0.047 ± 0.027 

0.038 ± 0.003 


0 042^® '^^^ 

0.02 


Table III: DO lepton asymmetries [50] 


III. FLAVOUR NON-UNIVERSAL AXIGLUON MODEL 

An axigluon is a massive, coloured ( SU{3)c adj), vector boson. Models of axigluon 
which have only axial couplings with the quarks have been suggested in the literature in 
many GUT like theories as chiral extensions of the QCD [57, 58]. Contribution to AFB 
for such a particle was studied even before the possible anomalous AFB was observed 
at the Tevatron in 2008 [10]. For this flavour universal, axially interacting massive 
gluon with coupling gs, the top quark AFB becomes negative for masses above niA ^ 
500 GeV. Upon the observation of a positive AFB by Tevatron in 2008, this model 
was found to be incompatible in the mass parameter regions allowed by the di-jet 
constraints from Tevatron. The AFB turns back positive if the assumption of universality 
of the interaction of axigluon with the quark families is dropped [12]. In our study 
here, we have used a more general, flavour non-universal axigluon with axial vector -i- 
vector couplings [13]. This model is obtained by breaking a larger symmetry group of 
SU{3)a X SU{3)b to the QCD colour group SU{3)c and a SU{3)c'- The axial-vector 
coupling of the axigluon to the first and second generation quarks is negative of that 
for the third generation and the vector couplings are the same for all three generations. 
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The couplings of the axigluon with quarks are described by the Lagrangian 

C = fyi‘T‘(^v+SA'r^)fA1„ ( 4 ) 

where T® are the Gell-Mann matrices. The couplings are parametrized by gv = 
" tanfze^) ’ Sa = sinfiOA) ’ generation of quarks. The parameters in this 

model are 9 a and We vary the value of the coupling in the range 6 a e [0, f ] which 
corresponds to varying the axial and vector couplings from a large value at small 6a to 
gv = 0,gA = gs for 0^ = f . A mass range of g [1,3] TeV is scanned. 

The decay width of axigluon and the density matrices for top-pair production me¬ 
diated by an axigluon are given in Appendix B1 . For an s-channel resonance, the 
terms in the tt pair production amplitude which are proportional to the linear power of 
COS0 (where 6 is the top-quark polar angle) contribute towards the AFB. The helicity 
dependent analysis of the top-quark decay distributions can give additional informa¬ 
tion about the bSM couplings. We will show in this study that this information can be 
accessed at the experiments from correlations among top polarization, top-quark and 
decay-particle asymmetries. 

We first discuss constraints coming from tt production cross section measurements, 
and top-quark level forward-backward and charge asymmetries measured at the Teva- 
tron and the LFIC (as appropriate). 

A. Constraints on the axigluon model 

We calculate the differential cross-section of the process {pp)pp tt ^ Ivht at the 
Tevatron with ^/s = 1.96 TeV and at the LFIC with = 7 TeV and a/s = 13 TeV for 
the SM -I- bSM with CTEQ6I [59] parton distribution functions with factorization scale 
fixed at Q = 2mt = 345 GeV, the top quark mass is taken to be mt = 172.5 TeV and 

as (mt) = 0.108. 

The cross-section calculated for the Tevatron, LFIC and the AC and AFB of the 
tt at those experiments in the axigluon model are shown in figure 1 . We constrain 
the model parameter space by limiting the predicted observables ApB 

and Ac to within la of the experimental values. As the values of 6a and niA grow 
larger, the couplings reduce, the mass of the mediating particle rises and the bSM 
contributions to the observables reduce. At large values of 6a, the figures correspond 
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Oa 



0A 


(a) Cross-section at the Tevatron with -y/s = 1.96 (b) AFB at the Tevatron with -y/s = 1.96 TeV 


TeV 




(c) Cross-section at the LHC with = 7 TeV (d) Aq at the LHC with y/s = 7 TeV 


Figure 1: Observables at the top quark level at the Tevatron and the LHC as a 
function of Oa for various values of masses for the axigluon. The experimentally 
measured values are marked in grey and the respective 2cr errors in dotted black 
lines. As the lines go from solid to dashed with larger gaps, the mass of the axigluon 
rises from 1 TeV to 2.7 TeV. 


to an axigluon model with only an axial coupling with the top quark and no resulting 
top polarization. For a lower mass range, constraints from the LHC allow only larger 
Oa and hence smaller coupling values, at the same time, interference with SM gives 
a constraint at the Tevatron which allows some region in the large coupling range as 
well. Ac gives a complimentary constraint and rules out large values of 6a (couplings 
close to gs) for a smaller mass of the axigluon. The result is that for the low masses 
of the axigluon, a range of couplings corresponding to 6a ^ (25° - 35°) and masses 
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niA ^ (1300 - 1900) GeV are allowed. Masses above these values are allowed for 
almost all parameter space with the only constraints coming from the Tevatron cross- 
section. 

CMS results constrain the mass of an additional massive spin -1 colour octet of par¬ 
ticles (eg. Kaluza Klein-gluon) which couple to gluons and quarks to above 3.5 TeV, 
which excludes the parameter region favoured by the experimental results from the tt 
process mentioned above [60]. The constraints can be evaded if the assumption of 
equal couplings of axigluons to light quarks and the top quark is relaxed. In this case, 
the values of coupling gv,gA we use can be split into gY,g\ and gv,g^ji^ where, the 
couplings with quarks would be constrained strongly from the axigluon direct produc¬ 
tion bounds. In the limit that the vector and axial couplings are equal or any one of 
the vector or axial coupling is small, our results can be recast into the modified model 
by using g^^^ = glf^glja- A more generalized version of such an axigluon model has 
already been discussed in the literature [61] along with constraints on the model from 
lepton and top quark asymmetries at the Tevatron and the LHC. 

The axigluon model can be constrained from B physics [62] results however, given 
the somewhat large hadronic uncertainties in some of the variables along with the 
possibility of relaxing these constraints in various modified axigluon models and/or 
by constructing UV completions, for the purpose of this study, we do not take these 
constraints into account. 

IV. U-CHANNEL SCALAR EXCHANGE MODEL 

In a second class of bSM models, AFB is explained due to contributions of a t- 
or M-channel exchange of new particles between the top quark-antiquark pair. The 
corresponding mediators do not show resonance behaviour and are elusive in the 
bump-hunting type analyses in tt pair production though they do contribute significantly 
to the angular distributions. We consider here a scalar particle called diquark, which, 
similar to a squark with R parity violation, transforms as a triplet under SU{3)c and has 
a charge of - 5 . The corresponding coupling is given by the lagrangian below, 

jC = t‘=T^{ys + ypy^)u(pa + h.c. (5) 

f = -ij^t* 
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( 6 ) 


= -KhVf 

We assume a right-handed coupling of the scalar with the up type quarks with y = 
ys = i/p. This ensures that flavour constraints and proton stability bounds are avoided. 
The density matrices for top pair production in the diquark model are given in Appendix 
B2. All calculations in this work are performed at tree level. The NLO contributions 
become important to study the effects on invariant mass distributions which we have 
not included here. These calculations are under progress for both axigluon and diquark 
models. 

A. Constraints on the coloured scalar model 

As in the case of an s-channel resonance in the previous section, the constraints 
are obtained from the measurements of the top pair production cross-section at the 
Tevatron and the LHC (7 TeV), the AFB and AC. We explore a parameter space of 
G [100,3000] GeV and y G [0, 27r] chosen so as to explore all the values of the 
coupling within the perturbative limit. As the value of the coupling rises, contribution of 
the bSM to all the observables becomes larger. For lighter diquarks, negative values 
of AFB and AC are predicted for large values of the coupling, though this mass range 
is ruled out by independent constraints from diquark pair production [63]. In figure 2 
we can notice from the top left panel that, as in the case of the axigluon, the Teva¬ 
tron cross-section provides the constraints in the parameter space of lower masses 
and couplings. In the next panel, the AFB measured at the Tevatron disallows lighter 
scalars and also constraints a part of coupling values for larger masses. The LHC 
cross-section constraints large coupling regions which give larger contribution and the 
cut-off coupling increases as mass of the scalar becomes heavier. The AC also allows 
larger coupling parameter space for higher masses of the scalar. 

The constraints from pair production of the coloured scalar from gluon fusion at the 
LHC are weak (-300 GeV) as reinterpreted from corresponding constraints on squarks 

[64] . There are further constraints on lower mass scalars from atomic parity violation 

[65] . Constraints from uu tt can be avoided by adding flavour symmetries (see for 
example [R"]). 
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TeV 



Yp (=ys ) 


(c) Cross-section at the LHC with = 7 TeV 



Yp (=Ys ) 


(b) AFB at the Tevatron with Vs = 1.96 TeV 



(d) Ac at the LHC with Vs = 7 TeV 


Figure 2: Observables at the top-quark level at the Tevatron and the LHC (7 TeV) as 
a function of the Yukawa coupling for various values of the diquark masses. The 
experimentally measured values are marked in grey and the respective 2(r errors in 
dotted black lines. The line spacing changes from solid to a dashed line with wider 
spaces as mass values rise from 100-2600 GeV. 


V. CONSTRAINTS FROM DUET PRODUCTION AT LHC 

The coloured scalar and vector bSM models get constrained from searches for di¬ 
rect production of bSM particles and subsequent decay to di-jet and four jet final states 
[qq ^ A^l] ,gg ^ (p(p^ 4;). Earlier constraints on axigluon model were obtained 

from searches of narrow resonances from dijet spectrum at 8 TeV LHC and were ex- 
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(a) Constraints on axigluon parameter space 



(b) Constraints on diquark parameter space 



Figure 3: Allowed parameter space for axigluon and diquark models are depicted as 
Green (lighter) coloured regions. Figure 3a shows the constraints from dijet searches 
as the blue (darker) shaded area. In figure 3b the dotted line represents the bound 
from pair production of bSM particles at LFIC. 


tended to the case of a wider width axigluon model with ^ = 0.3 where the axigluon 
has only axial-vector couplings [67] (also see [c-j for bSM particle off-shell effects in 
dijet searches). We reinterpret these constraints to the case of the axigluon model in 
this study which has both axial-vector and vector couplings with the quarks and find 
the excluded parameter range where ^ < 0.3. Figure 3a shows the parameter space 
allowed for the axigluon model. The following constraints are put on the model param¬ 
eters to obtain the allowed values : reinterpreted searches for bSM resonances in dijet 
production, tt cross-section and top charge asymmetry measurements at 7 TeV LFIC 
and cross-section, top forward backward asymmetry measurements at the Tevatron as 
discussed in section III A. The coupling values corresponding to 9a > 27° are ruled out 
for axigluon masses up to ^ 4 TeVs as these narrow, resonant particles would have 
been detected in the dijet searches. The allowed values of couplings correspond to 
6a ^ 10° - 27° for the mass range between 1.5 TeV to 3 TeV. Note that the constraints 
from dijet searches may be relaxed if the magnitude of coupling of axigluon is different 
for the third generation of quarks as compared to the first and the second generations. 
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The diquark mass is bound from below to 300 GeV from pair production of 

diquarks via gluon fusion at the LHC [64]. The direct production bounds along with 
the constraints obtained from top quark pair production cross-section and top charge, 
forward backward asymmetry measurements at LHC and Tevatron (see section IV) 
are shown in figure 3b. A narrow strip of parameter space is allowed when couplings 
are large due to destructive interference effects. Besides this region, the rest of the 
allowed diquark parameter space follows the expected behaviour of small coupling 
values < 0.2 for lower masses and for a diquark of mass 3 TeV, xjs as large as 2 is 
allowed. 

VI. POLARIZATION OF THE TOP QUARK AND DECAY-LEPTON DISTRIBUTIONS 

The decay kinematics of leptons embeds the information regarding top quark pro¬ 
duction dynamics, kinematics and polarization [20]. Different lepton observables em¬ 
bed these effects in different ways and so provide a number of probes which are all 
correlated with the top quark kinematics and polarization. For a detailed analysis of top 
quark decay see [31 , 32]. In this section we discuss distributions of the lepton polar an¬ 
gle, azimuthal angle and energy in SM decay of top quark and construct asymmetries 
based on these distributions to probe top quark bSM interactions. 

A proper treatment of the decay distributions of the top quark requires the spin 
density matrix formulation, which preserves correlations between the spin states in 
the production and in the decay. 

The spin density matrix for t in the production of a tt pair with the spin of t summed 
over, can be expressed as 

P tt production production {^t/^t') production C^) 

Af 

The density matrix gets SM contributions ^“^cl respectively 

from gluon-gluon and quark-anti-quark initial states, a contribution 
the bSM model, and a contribution pinterference {^t,K) from the interference between 
the SM amplitude and the bSM amplitude: 

P {At, + Pgj^ (Tf/Af) + PJjSM {^t/Aj{^ + pinterference {^t/Aj.^ (8) 
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The spin density matrix for the decay of the top quark is given by 


^top decay decay 


■^decay {K) > 


( 9 ) 


with the spins of the decay products summed over. 

The squared amplitude for the combined process of production and decay is given 


by 


M\^ 


m}) 

Ttmt 


Ep(Af,A;)r(Af,A;) 

Af,A' 


( 10 ) 


This expression assumes a narrow-width approximation for the top quark. Top decay 
is assumed to progress through SM processes. In the rest frame of top quark, the 
differential decay distribution of the top quark is given by 

dTt 1 + ApkfCosO 

-^— =- — - ( 11 ) 

Tdcos{6) 2 


where 6 is the angle between top-quark spin direction and the momentum of the decay 
product /. For N(At) number of top quarks with helicity Aj, polarization Ap is defined 


N(Aj = +)-N(Af = -) 
N(Af = +) + N(Af = -)' 


( 12 ) 


and the coefficient kf is called the top-spin analysing power of the decay particle /. 
For the case of leptons as the final state particles, the factor kf = 1 at tree level in the 
SM. When the top quark is boosted in the direction of its spin quantization axis, eqn 


(11) gets modified to 


drBoosted _ (1 - (1 + Af cos Oil - j6(cos Oil + At)) 

rdcos(0t/) 2(1 - |6cos0f;)3 

where Oa is defined as the angle between lepton and top quark momenta in the 
boosted frame. Lepton kinematic distributions for the tree-level SM differential cross- 
section for the process pp^tt^l+ jets with ^/s = 1.96 TeV are presented in the 
figure 4. The energy and azimuthal lepton distributions are uncorrelated with the po¬ 
larization in the rest frame of the top quark, though correlate with the polarization in 
the boosted frame. Fligher-order corrections to the production and decay processes 
have been calculated for the SM and the distributions are found to be qualitatively 
unchanged [69, 70]. Due to this reason, we expect that the effect of higher-order 
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(a) Lepton polar angle distribution In top quark rest (b) Lepton polar angle distribution In the lab frame, 
frame. 




(c) Lepton Azimuthal angle distribution In the lab (d) Lepton energy distribution In the lab frame 
frame 

Figure 4: The tree-level lepton polar and azimuthal distributions for 
pp^tt^l+ jets with ^/s = 1.96 TeV. In the above plots, the average boost of the tt 

pairs is 0.34. 

corrections to the asymmetries constructed from decay-lepton distributions to be rela¬ 
tively small as the corrections partially cancel out within the difference and ratio taken 
to derive the asymmetry. 

A lepton polar-angle asymmetry with respect to the top-quark direction can be de¬ 
fined by 

Ati _ ^(cos(gtf) > 0) - (r{cos{9ti) < 0) 

o-{cos{eti) > 0 ) + o-{cos{eti) < 0 ) 

has been measured at both LHC and Tevatron in the lab and tt center of momen- 
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turn frame, albeit with large statistical errors and different results from CDF and DO 
[49]. Integrating eqn (11) the top rest-frame lepton asymmetry can be related to the 
polarization of the top quark, 


= (15) 

In QCD, A%^^^ = Ap = 0, though in the boosted frame, the lepton polar asymmetry 
with respect to the top quark is large even in a tree-level SM calculation. 

In the lab frame where the top quarks and leptons are boosted and the cross- 
section convoluted with the pdf, the correlations between various angles and energies 
become more complicated. The lepton polar angle with respect to the proton beam 
is a convenient observable which does not require top-quark rest frame or momenta 
reconstruction. The lepton polar asymmetry A^g in the lab frame is also 0 at tree 
level in SM QCD. Aj,g is identically 0 at the LHC due to the symmetric nature of the 
initial state. This asymmetry according to our analysis correlates the best with the off- 
diagonal elements of the top quark density matrix for the I + jets process considered 
(see Section VIII) for both axigluon and diquark models. An analytic study of the lepton 
polar angle and its correlation with top AFB and polarization has also been made by 
Berger et al. [71]. They relate the lepton and the top quark level polar asymmetries 
and subsequently use this relation to distinguish between a sequential axigluon and a 
W’ type model [72]. 

In the top-quark rest frame, other lepton kinematic variables : azimuthal angle and 
its energy have no dependence on the helicity of the top quark and hence the inte¬ 
grated asymmetries are uncorrelated with the polarization. It has been noted in the 
literature that the lepton azimuthal distributions correlate with the polarization of the 
top quark in a boosted frame [32, 73-75]. Sums and differences of azimuthal decay 
angle in top pair production process have also been used in the literature to study the 
polarization and spin correlations of top quark in detail [76]. For a detailed analysis of 
analytic relation between polarization of a heavy particle and decay particle azimuthal 
asymmetry, see [77]. We reproduce the azimuthal distribution in the lab frame for the 
SM tt pair production process at the Tevatron in figure 4c. The azimuthal distribution 
can be measured at both Tevatron and LHC and requires only partial reconstruction of 
top quark rest frame. The azimuthal angle is defined by assuming that the top quark 
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lies in the x-z plane with proton (beam) direction as z-axis. From this distribution, an 
azimuthal asymmetry about a point (po can be defined as 



(16) 


A natural choice for the value of (po would be the point of intersection of the distributions 
corresponding to left and right helicity top quarks. For SM, this point is about (p = 50° 
for both Tevatron and the LFIC. The SM point would correspond to 0 polarization and 
would maximize correlation with bSM contribution. We assume a value of (po a bit lower 
at 40°. Since the positive helicity top quark have larger differential cross-section in 
this region, this choice enhances correlations of the lepton level asymmetry for larger 
positive (or smaller negative) values of polarization. The standard model tree-level 
results for this asymmetry at the Tevatron and the LFIC respectively are given in table 
IV . In the lab frame, due to the boost and rotation from the direction of the top quark, 
Acf, is sensitive to both the polarization and the parity breaking or t-channel structure 
of the top quark coupling. Another observable which can be constructed from the 
decay-lepton kinematics is the lepton energy asymmetry about a chosen energy Eo'. 


^ tr(E; > Eq) - (r{Ei < Ep) 

cr(E/ > Eq) + cr{Ei < Ep) 


(17) 


No reconstruction of the top-quark rest frame is needed to measure Ei. Just like the 
azimuthal case, this asymmetry can be measured both at the LFIC and the Tevatron. 
The lepton energy distribution is sensitive to the polarization of top quark [32], as 
shown in figure 4d. Similar asymmetries based on the energy of decay particles or the 
ratios of these energies have been used in the literature to study bSM physics [74, 78- 
80]. We define the lepton energy asymmetry about a value of Ep = 80 GeV, to act 
as a better discriminator between bSM and SM. Ideally, the point of intersection of the 
positive and the negative top-polarization curves should form the best correlation with 
the top polarization, though this point varies with the energy and the invariant mass of 
the initial state. Standard model values of asymmetries mentioned in this section are 
given in table IV. It would be interesting to use SM distributions at NLO to decide the 
reference points Eo,(po, but since in the end we construct asymmetries, we expect that 
the qualitative behaviour of our results would not change. 

In the recent past, polarization measurements have been made by collaborations 
both at the Tevatron and the LFIC. The polarization at the Tevatron points towards a 


18 




Asymmetry 

Q = mt 
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^FB 

0.645 

0.642 


-0.113 

-0.116 

^E, 

0.381 

0.397 


(a) Tevatron a/s = 1.96 TeV 


Asymmetry 

Q = mt 

S 

(N 

II 

o 

Atl 

^FB 

0.748 

0.748 


-0.075 

-0.077 

^E, 

0.138 

0.146 


(b) LHC = 7 TeV 


Asymmetry 

Q = mt 

S 

(N 

II 

o 

Atl 

^FB 

0.789 

0.788 


-0.041 

-0.044 


0.036 

0.038 


(c) LHC Vs = 13 TeV 

Table IV: Scale dependence of SM values of various asymmetries tree level. 


small positive value and that at the LHC to small negative values. This is consistent 
with the small coupling and large mass regimes of both the models studied here. 

In the next section, we use correlations among top charge and forward-backward 
asymmetries, decay lepton angular and energy asymmetries, and polarization to un¬ 
cover specific properties of bSM particles which can be inferred from the Tevatron and 
the LHC data. 

VII. CORRELATIONS 

The parameter space of c [1000,3000] GeV and 6a ^ [10,45] are explored for 
the axigluon model and c [100,3000] GeV and i/s e [0,27r] for the coloured scalar. 
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The figures in the section VIIB show parameter space allowed by the constraints men¬ 
tioned in sections III, IV. 


A. Correlations between charge and forward-backward asymmetries 


The correlation between the AC at 7 TeV LHC and the AFB at the Tevatron have 
been used in the literature constrain various bSM models(see for example [14]). These 
constraints are model dependent and the asymmetries are not in general tightly cor¬ 
related [81]. We show similar correlations in figure 5 where we plot Ac vs Afs, using 
the relation 

A NLO I AbSM M 

^C/FB-^C/FB +^C/FB 

This relation is valid as long as the bSM physics corrects the SM cross-section of the 
tt pair production process by a small amount. 



0.04 

0.02 

^ 0.00 
- 0.02 
- 0.04 


0.00 0.05 0.10 0.15 0.20 0.25 0.30 



A 


t 

FB 


(a) The red markers represent axigluon model with (b) diquark model with right handed couplings to 
only axial interactionsfg^y = 0,^^ = gs). The size u,t quarks. The size of the plus marks represent a 
of the plus marks represent a mass range from mass range from 100-3000 GeV 


1000-3000 GeV 


Figure 5: Correlation between top-quark asymmetries vs at the Tevatron and 
the LHC (v^s = 7TeV). The grey solid and dashed lines represent the observed 
values of the respective asymmetries and their la errors. 


The pure axial axigluon which leads to unpolarized top quark is disfavoured as it 
does not have a parameter space where it can explain both Afg and Ac experimental 
values. In the diquark model, the coupling to right handed quarks is sampled from 0 
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(a) Lepton polar asymmetry Tevatron = 1.96 (b) Lepton polar asymmetry Tevatron Vs = 1.96 

TeV. TeV. 

Figure 6: Correlations between top AFB and lab frame 9i,9ti asymmetries at the 

Tevatron ^/s = 1.96 TeV 

to 271 where large mass or small couplings lead to a better agreement with SM NLO 
values of asymmetries. 

B. Correlations among lepton and top asymmetries 

In this section, we study the correlations among top polarization, top asymmetries 
and decay lepton asymmetries. We show that combined, they form sensitive dis¬ 
criminators between models with different dynamics. The top-quark and decay-lepton 
asymmetries are calculated at various points in the parameter space allowed by the 
experimental constraints discussed in section III and IV. The expected polarization 
of the top quark, for corresponding points in the parameter space, is represented in 
colour contrast form inside the graphs and clear trends for the polarization can be ob¬ 
served. In all the following figures, top-quark asymmetries represented on the x-axis 
are calculated as shown in eqn (18) and the lepton asymmetries shown on the y-axis 
are calculated including SM-i-bSM contributions at tree level. 

1. Asymmetry correlations for the Tevatron 

The correlations of the lepton-level asymmetries with the top AFB at the Tevatron 
are shown in figures 6a, 6b, 7a and 7b. For the case of axigluon as its mass is in- 
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(a) Lepton Azimuthal asymmetry about (po for the (b) Lepton energy asymmetry about energy Eq for 
Tevatron yi = 1.96 TeV. the Tevatron -y/s = 1.96 TeV. 


Figure 7: Correlations between top AFB and the lepton energy and azimuthal 
asymmetries at the Tevatron ^/s = 1.96 TeV 

creased, the polarization rises until ^ 1650 GeV and then drops again for even 
larger masses. The diquark model predicts negative polarization for a significant por¬ 
tion of parameter space, turning positive only for large couplings. The large mass 
region for the diquark also favours large (negative) values of azimuthal asymmetry, 
smaller lepton polar asymmetry and larger lepton energy asymmetry. Lepton polar 
asymmetry correlation with top AFB shows large overlap between the two models. 
The observed value for the lab-frame lepton polar asymmetry in tables 11,111 points to¬ 
wards a positive polarization between = 0.2 to 1.1. In this region (see figure 6a), a 
large positive value for polarization is favoured for the diquark model and a large con¬ 
tribution to top AFB. The axigluon model is compatible with both the observed value 
of Ae^ and a small contribution towards longitudinal polarization for a significant part 
of its parameter space. Figure 6b shows the asymmetry in the lepton polar angle with 
respect the top direction, Ag^^, which is equal to twice the top polarization, eqn (15), 
when calculated in top-quark rest frame. It receives contribution from bSM physics via 
the boost of the parent top quark. In the lab frame, large deviations of Ag^^ from the 
SM value correlate with large contribution to the top-quark polarization from the bSM. 
The values of asymmetries grow closer to the corresponding SM values with increase 
in mass and reduction in the bSM coupling strength. The azimuthal asymmetry and 
lepton polar asymmetry with respect to the top-momentum direction {Ag^^) and lepton 
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(a) Lepton Azimuthal asymmetry about = 40 
degrees for the LHC ^/s = 7 TeV 
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(b) Lepton polar asymmetry(w.r.t top quark 
direction) for the LHC Vs = 7 TeV 



(c) Lepton Energy asymmetry about Eq = 80 GeV 
for the LHC Vs = 7 TeV 


Figure 8; Correlations between lepton and top kinematic asymmetries at the LHC-7 

TeV. 


energy asymmetry in figures 6b, 7a and 7b discriminate well between the s-channel 
and u-channel exchange models though the parameter spaces within the model are 
clumped together. When combined with polarization, all correlations enhance their dis¬ 
criminating power especially to distinguish between s-channel and w-channel models 
as they predict opposite signs of polarization for a large portion of parameter space. 


2. Asymmetry correlations for the LHC 

The lepton-level asymmetry correlations with tt charge asymmetry are shown in 
figures 8a-8c for the LHC 7 TeV run and figures 9a-9c for the LHC 13 TeV. The plots 
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O-Axigluon D-Diquark 



Figure 9: Correlations between lepton and top kinematic asymmetries at the LHC-13 

TeV. 


are made for the region of the model parameter space constrained in sections III and 
IV. For the = 7 TeV calculation, we use mt = 172.5 GeV and factorization scale 
Q = 2mt and rtg = 0.108 to remain consistent with the ATLAS and CMS reconstruction 
of Ac- For ^ = 13 TeV,, the mass of top quark is chosen at the updated central value 
mt = 173.2GeV and Q = 2mt with ccs = 0.108 and CTEQ6I pdf. 

A large portion of the parameter space predicts a negative polarization at 7 TeV LHC 
for the axigluon model. The diquark model predicts a small negative polarization for 
small couplings with quarks. When heavier diquark models are considered, larger cou¬ 
plings are allowed leading to a large positive contribution to the polarization. Observed 
values of polarization from CMS and ATLAS are compatible with -0.03 < Ap < 0.07, 
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which covers a large region of parameter space for both axigluon and diquark models. 
As in the case of the Tevatron, polarization is an important discriminant between mod¬ 
els for the LHC as well, especially when combined with decay-lepton asymmetries. It 
is able to distinguish overlapping parameter space regions between the two models. 
This is more true when the couplings are small and the bSM effects are more difficult to 
detect as the s-channel and u-channel exchanges predict small polarization, but with 
opposite signs in this region. The energy asymmetry becomes smaller for the LHC at 
13 TeV due to the effect of the overall boost. The values of azimuthal and polar asym¬ 
metries do not change significantly for higher energy and so remain good observables 
for the study of top quark dynamics. 


VIII. ASYMMETRY CORRELATIONS AND TOP TRANSVERSE POLARIZATION 


As remarked earlier, keeping full spin correlations between the production and de¬ 
cay of the top quark in a coherent manner requires the spin density matrix formalism. 
In this formalism, the top polarization, which played a significant role in the above 
analysis, corresponds to the difference in the diagonal elements of the density matrix, 
as seen from eqn (12). The off-diagonal elements of the density matrix can also be 
significant in practice, and they would contribute to the transverse polarization of the 
top quark, corresponding to a spin quantization axis transverse to the momentum. In 
SM, these terms arise at loop level and have been studied in the literature along-with 
transverse polarization and observables have been suggested to measure their con¬ 
tribution [82]. We examine in this section what role these off-diagonal matrix elements 
and transverse top polarization play in the two models considered in this study. 

Following the formalism developed in [32], the spin density matrix integrated over a 
suitable final-state phase space can be written as ( 7 (A,A') = where crtot 

represents the unpolarized cross section. The matrix Pt(A,A') can be written as 




I 1 + 73 71 - ^’72 
y 7i + ^’72 1-73 


(19) 


Here 73 is the longitudinal polarization, j/i and j /2 are polarizations along two transverse 
directions. The expressions for the//; in terms of the top-quark density matrixes (At, A+ 
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can be written as, 


V3 = 

m = 

if]2 = 


[cr 




O'tot 

{cr{+-) + a{-+)) 
O'tot 

((r(+-) -C7(-+)) 
O'tot 


( 20 ) 

( 21 ) 

( 22 ) 


Splitting the top density matrix as shown in eqn (10) under the narrow-width approx¬ 
imation, the helicity-dependent decay density matrix in the rest frame of top quark 
separates into a simple functions of the decay angle: 


dr(A,A') = c X A{\,\')dO,i 


(23) 


where 


a;. A' ^ 

+ 

+ 

l + cos(0;) 

— 

sin(0/)e“*^' 1 —cos(0/). 


O/ is the solid angle in which the lepton is emitted and c is the integrated contribution 
of the rest of the decay kinematic variables. The resulting lepton angular distribution 
in the lab frame is , 

dcos^l^i) d(pi ^ ^ ^ ^ 

The off-diagonal elements in the top-quark production density matrix do not contribute 
to the total cross-section due to an overall factor of sin (0/) which integrates to 0. They 
do contribute instead to the kinematic distributions of the decay particle, although this 
effect is quite small for most observables. 

In this study, we find that the lepton polar angle asymmetry defined in the lab frame 
is sensitive to the off-diagonal terms in the top quark density matrix eqn (7). The 
transverse polarization originating from these off-diagonal terms contains further infor¬ 
mation about the dynamics of top-quark interaction. This relation has been pointed out 
before in the context of a wide-width colour octet bSM particle [83, 84]. 

In figure 10 we study the contribution of off-diagonal terms to the lepton distributions 
and present the distributions for a few sample masses of bSM particles. 

It can be seen that the contribution of the off-diagonal density matrix elements can 
be significant, and is particularly important for the diquark model. These can in turn 
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(a) axigluon 


(b) diquark 


Figure 10: Contribution to Lepton asymmetry from the off-diagonal terms of top 
density matrix calculated in the lab frame for the Tevatron ^/s = 1.96 TeV for the 
axigluon and diquark models. The red lines represent the asymmetry for diagonal 
density matrix and the black line represents the distribution for the case of total 
density matrix. The darker lines represent allowed regions of parameter space. 


lead to significant transverse polarization of the top for appropriate range of parameters 
which could be measured experimentally. 


IX. CONCLUSIONS 

The forward backward asymmetry of the top quark in top-pair production process 
at the Tevatron collider was, for a long time, anomalously large and a persistent effect 
observed independently by both DO and CDF detectors. It has been demonstrated 
only relatively recently that NNLO contributions give rise to an of the right order 
of magnitude and seems to be in agreement with the values measured experimentally. 
Previously, many bSM models had been proposed with parity breaking interactions to 
explain the observed ^FB- Many of these models predicted a charge asymmetry at 
the LFIC. Since LFIC has a gluon dominated initial state as opposed to the Tevatron 
where qq was the primary initial state, the asymmetries predicted for LFIC coming from 
the bSM couplings to the quarks get diluted. The data gives values for Ac consistent 
with the SM and so far there has been no evidence for the new particles predicted 
in the different bSM models. Linder these circumstances there is a need to construct 
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measures which can distinguish between different sources of the either SM or 
bSM. 

One such measure is provided by polarization of the top quark which has a non¬ 
zero value in the presence of a parity breaking interaction. Within SM, top polarization 
is close to 0. Observables that correlate with top polarization can be used to distin¬ 
guish between various SM and bSM contributions. In continuation to a previous work 
where correlations between polarization and forward backward asymmetry were used 
to constraint bSM [14], we have introduced the correlations between lepton polar, az¬ 
imuthal and energy asymmetries and top charge asymmetry and showed how they 
can be used together with top longitudinal polarization to distinguish between SM and 
bSM. 

In the reference [85] the authors have constructed dilepton central charge and az¬ 
imuthal asymmetries and studied it along with top quark polarization, forward backward 
asymmetry and tt spin correlations for benchmark models of G’ and W’. Subsequently, 
in reference [71] the authors have shown that the lepton polar asymmetry and top for¬ 
ward backward asymmetry and lepton charge asymmetry vs top charge asymmetry 
correlations can be useful in the study of W’ and G’ models. Our work adds multiple 
new observables to the analysis of new physics in tt pair production which include 
single-lepton azimuthal angle, energy and polar angle (wrt top quark) in the lab frame 
which show signatures from parity breaking in top interactions and help isolate and 
constrain the interactions of bSM particles. 

We demonstrate the efficacy of the correlations between forward backward asym¬ 
metry and the lepton asymmetries at Tevatron and charge asymmetry and lepton an- 
gluar,energy asymmetries at LHC, by utilizing a representative s-channel model, ax- 
igluon and an u-channel model, diquark. Constraints on these models are obtained 
based on measured values of tt cross-section at Tevatron and LHC 7TeV, and 
A^ and resonance searches in dijet, four jet cross sections. The parameter space 
of axigluon allowed within 2-o- of the measured values of the stated observables in¬ 
cludes a lower bound on the mass of axigluon at 1.5 TeV with a corresponding cou¬ 
pling 9a > 27°. The allowed mass of the diquark is bounded from below by 300 GeV 
and masses above are allowed with the coupling of the model bounded from above by 
a value of 0.2 for smaller masses which rises to i/s < 2 corresponding to = 3 TeV. 
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Another sliver of parameter space is allowed for larger couplings of the diquark due to 
destructive interference effects. 

For the first time we have presented the complete density matrix of the top quark, 
including the off-diagonal elements for top quark pair production process, in the ax- 
igluon and diquark models to aid further studies. We use these to show that the lepton 
polar asymmetry in the lab frame shows a correlation with the transverse polarization 
of the top quark for axigluon model and even more significantly for the diquark model. 
The lepton asymmetry usually considered in studies of top polarization is calculated 
from lepton polar angle with respect to the top quark and it does not show this correla¬ 
tion with transverse polarization. The correlation of transverse polarization with lepton 
azimuthal or energy asymmetry is also very small. 

Finally, we extend our analysis to 13 TeV LFIC where, even though the values of the 
asymmetries get diluted, the correlations between accurate measurements of charge 
asymmetry and lepton asymmetries still separate out bSM and the SM in the 2 di¬ 
mensional space. Taken together, these correlations can indeed be used to improve 
significance of the constraints on bSM from LFIC data even in the initial stages of low 
luminosity. 
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Appendix A: AC at 13 TeV LHC 

Since the AC calculated at NLO for 13TeV LFIC was unavailable at the time of 
submission of this work, we note that the available charge asymmetry values [48] form 
a smooth function of the beam energy and fit them to a polynomial to find the AC as 
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a function of the beam energy. We obtain a fit to a polynomial presented in eqn (A1) 
with goodness of fit parameter = 0.9995. 


Ac (^/s) = 3.12 X 10“2 - 4.37 x 10“^+ 2.6269 x 10“^s - 5.683 x 10“S2 (A1) 

This gives a value of Ac(13 TeV)=0.0063. 


Appendix B: tt Production density matrices 


1. axigiuon density matrices 


With Q = cos(0t), So = sin(0t), = -y/l - ^ and = W1 - 
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X {^gA {gv - gA^) Ce + gv - 4 ) + ^ (^A - gv^) C 20 ) ) 

(B8) 

To present the dependence on top boost and polar angle clearly the amplitude 
square is written in terms of the polar angle 9 in tt center of momentum frame. The 
off-diagonal terms in the gluon initiated process are zero and the diagonal terms in 
gluon initiated process are not dependent on the top quark polarization therefore we 
have omitted these here and they can be found in many references including [14]. 
Decay width of axigluon at tree level is given by, 

Ta = ^^{gA {^A i^A + 5) - + gl (m\ {^a + 5) + } (B9) 

2. diquark density matrices 

The top-quark spin density matrix for the w-channel exchange is given below in tt 
center of momentum frame. The notation and SM only contributions remain the same 
as for the case of the s-channel model. 
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